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Abstract 

We compare several recent theoretical studies of the single transverse spin asymmetry in dijet- 
correlations at hadron colliders. We show that the results of these studies are all consistent. 
To establish this, we investigate in particular the two-gluon exchange contributions to the relevant 
initial and final state interactions in the context of a simplifying model. Overall, the results confirm 
that the dijet imbalance obeys at best a non-standard or "generalized" transverse-momentum- 
dependent factorization. 
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1. Introduction. Recently, there has been tremendous interest in the "imbalance" of 
two jets that are produced nearly back-to-back in azimuthal angle in a hadronic reaction. 
Especially the single transverse spin asymmetry (SSA) for this process has received a lot 
of attention in theoretical work [2I, I^, 0, Q, Q, [s], ^, while experimental studies have 
begun at the Relativistic-Heavy-Ion-Collider (RHIC) at Brookhaven National Laboratory 
(BNL) Q. In Ref. it was proposed to use dijet-correlations at RHIC to learn about 
the transverse-spin and transverse momentum dependent (TMD) Sivers functions [11]. A 
transversely polarized proton (with momentum Pa and polarization vector S_[_) is scattered 
off an unpolarized proton with momentum Pb, producing two jets with momenta Pi and 

A{Pa, S^) + B{Pb) ^ Ji(Pi) + J2{P2) +X . (1) 

For convenience, two additional momentum vectors, P = (Pi + P2)/2 and q = Pi — P2, 
are introduced. The transverse momentum components of P and q are relevant for our 
discussions: Pj_ is the average of the transverse momenta of the two jets and q± represents 
their imbalance. The most interesting kinematic region arises when the imbalance between 
the two jets is small: q± ^ Pi, so that there are two very separate momentum scales in the 
problem. 

Much of the very recent theoretical work has centered on the question if cross sections 
'or the process may be factorized in terms of TMD functions in this kinematic regime 

mm 

. This question has been addressed by three groups, in the following ways: 



(1) inRefs. , the gauge-link structure of the TMD parton correlators appearing in the 
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,ll3|, 



14, hj that 



dijet process was investigated. It had been known beforehand 
Sivers-type SSAs in semi-inclusive deep inelastic scattering (SIDIS) or the Drell-Yan 
process may be generated by final-state or initial-state interactions, respectively, of the 
involved partons from the polarized hadron. These interactions sum to a simple gauge 
link in both cases, which carries the color charge of that parton, albeit - famously 
- with an opposite sign between the two processes. As was discussed in , the 
SSA for dijet (or dihadron) production is associated with both initial-state and final- 
state interactions of the partons involved in the hard-scattering. It was found that if all 
interactions are summed up, the resulting gauge link for the TMD parton distributions 
takes a much more complicated form. It does not involve just the color charge of the 
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relevant parton, but has in general knowledge about the full hard-scattering process 
and its color structure. As such, different correlators were found to appear in different 
partonic channels, even if the parton type entering from the polarized proton is the 
same. This observation makes the non-universality of the TMD parton distributions 
much more dramatic than previously indicated by their sign difference between the 
SIDIS and Drell-Yan processes. It was also found in I2I, that a certain weighted 
moment of the spin-dependent cross section reduces the expression to a more standard 
form akin to coUinear factorization, involving the twist-three matrix elements of 16, l3] 
and special hard-scattering factors. 

(2) in Refs. 0, H], a different approach was taken. The SSA for the dijet imbalance was 

examined at first order in perturbation theory, starting from collinear factorization 
in the intermediate transverse momentum region Aqcd <^ Qi. ^ P±i and carefully 
examining the one-gluon radiation contributions. It was found that, at this order, the 
cross sections do factorize into TMD parton distributions that follow their definitions 
in either the SIDIS or the Drell-Yan process, while all the initial-state and final-state 
interaction effects can be absorbed into the hard-scattering factors. These leading- 
order hard factors can be calculated either in a model-inspired fashion or in a partonic 
scattering picture \7\. They coincide with the ones found in for the weighted 

spin-dependent cross section mentioned above. 

(3) Ref. {s] stresses that the non-trivial gauge hnk structure found in 0, implies that a 

"standard" factorization in terms of TMD parton distributions cannot hold for this 
process. A particularly transparent example is given in terms of a simple abelian 
model in which the two scattering partons are assumed to have different and unrelated 
charges, gi and g2. By a first-order calculation it is suggested that indeed the gauge 
link to be associated with the parton of charge gi will in general depend on the charge 
g2 as well, in violation of standard factorization. From a phenomenological point of 
view, TMD correlators extracted from dijet- correlations would have no connection with 
those from the SIDIS and Drell-Yan processes, because their definitions are different. 
It is argued that if a standard definition of parton distributions is kept, there will be 
uncancelled singularities at higher order of perturbation theory. This happens for the 
spin-dependent as well as the spin-averaged cross section. 
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The goal of the present paper is to put the results of Refs. into context . 

It is clear from the above that there is not necessarily a contradiction among the results. 
(1) and (3) agree in their assessment that universality, and hence standard factorization, of 
the TMD distributions is broken for this process. (2) is based on a first-order calculation, 
and it is conceivable that deviations from the factorized structure found in Refs. 0, Q] 
become apparent only at higher orders. In order to shed more light on this, we perform the 
following studies: first, we compare the results in (1) and (2) by expanding the gauge links 
given m to first order in the strong coupling. This should lead to the results of (2), 

if there is mutual consistency. Furthermore, we will consider the two-gluon contributions 
to the initial and final-state interactions. This will unambiguously show if the gauge link 

nn n 

exponentials have the complicated and process-dependent form presented in [2|, |3| and [8| , 
which would demonstrate that indeed standard factorization breaks down beyond the order 
considered in (2). It will also address the consistency between (1) and (3). For siniplicity, 
we will perform the calculation in the context of the model considered in Ref. {2!, Q], and 
only for one underlying partonic process, qq' qq' . 

2. First-order expansion of the gauge-links in Refs. 

aa. 



Ref. 



3(] finds the following gauge link in the correlator for the partonic channel qq' —>■ qq' 



contributing to the SSA for the dijet imbalance: 



m - 1 



(2) 



which is normalized to the correlator for unpolarized scattering. Replacing the factor 
1/{N^ — 1) with 1/(AN^) will recover the normalization used in In the first-order 

perturbative expansion, we will find the following contributions from the various gauge links 
in the above equation, 

1 



-k+ + ie 
1 



1 



(3) 



-A;+ + ie A;+ + ie 

Here we have suppressed the color matrices. Since the latter are traceless, one has Tr (WP^) = 
1 to this order. Using these expansions, we get the following first-order expansion of the 



^ We will from now on use for simplicity the labels "(1),(2),(3)" to refer to the papers @,0|i and 

respectively, following the list given above. 



gauge link 



qq'- 



Uqq' 



^ - 3) ' 



k+ + ie 



+ (-2) 



k+ + ie 



Comparing with the results in 



(4) 



we fine t^iat the first term corresponds to the sum 



of the two final-state interaction factors of [a, l7|, and the second term to the initial-state 
interaction factor there. More precisely, the coefficient of {—ig)/i—k~^ + i^) coincides with 
the factor {Cpi + Cf2)/Cu found in column (1) of Table I of [7|, while the coefficient of 
{—ig)/{k^ + ie) is identical to Cj/Cu- This comparison can be extended to the diagrams for 



the other partonic channels, and agreement between the results of 



31] and 



is found 



{zg)/{k^ 



in each case. In addition to Eq. ([3]), one then also needs W'^ 
3. Two-gluon exchange contribution. 

In this section we study the two-gluon exchange contribution to the initial and final-state 
interaction effects in the dijet imbalance. As we stated earlier, we will for simplicity use the 
abelian model of 0, 8| in which the two scattermg partons are assumed to have different and 
unrelated charges, gi and g2- We also follow [8] to consider only the case of an underlying 



qg 



qq' hard process; however, we take the q and q' as fermions and not as scalars as in 



To set the stage, we write the differential cross sections for the dijet-correlation in this 
model in the factorized forms assumed in 

dAa{S^) e^'^S'lql 



dyidy2dPld'^q_i 



^1 



d ki±d k2±d X_ 



X X, qfriXa, h^) X, qr'^X,, k2^) 



Mr 



(5) 



X 



for the transverse-spin dependent case, and 



dyidy2dPld'^q± 



d:'ki±d:'k2±d:'X± Xaqf^^^ixa, ka±) Xtq^'^'^'^iXb, kb±) 



SIDIS / 



X [S,,,^,,iX^)H;^,^^^,{Pl)]^ 6^'\ka± + h± + X^- gl) , (6) 



for the spin-averaged case. Hqqi^qqi and Sqqi^qqi are partonic hard and soft factors, respec- 
tively, and the [ ]c represents a trace in color space between the hard and soft factors due 
to the color flow into the jets [l8|, [l9|. The hard factors in Eqs. ([5]), ([6]) only depend on the 
single hard scale in terms of partonic Mandelstam variables of the reaction qq' —>■ qq', 



and Xa = ^ (e^^ + e^^), = ^ (e~^i + e"^^) with yi and ?/2 the rapidities of the two jets. 
In Eq. ([S]), q^^^^ and g^iDis (jgj^ote the transverse-spin dependent Sivers quark distribution 
for hadron A and the unpolarized TMD quark distribution for hadron respectively; these 
TMD parton distributions were chosen to follow their definitions in the SIDIS process. 

In the present model where the two initial partons have separate charges, the associated 
parton distributions will depend on these charges. For example, for a polarized hadron 
with momentum P4 = (-P4,0~,0_l) with = l/v^(P° ± P4) and transverse spin vector 
S'_L, the TMD distribution for quark flavor a with charge gi can be defined through the 
decomposition of the following matrix element. 



(7) 

where the gauge hnk C is defined in a covariant gauge as 

^vaigu 00; ^) = exp (^-igi dXva ■ A{Xv + , (8) 

with the path extended to +00. Va is a vector conjugate to the momentum vector Pa- Since 
we will work in a covariant gauge throughout this paper, the vector Va could be chosen to 
be a light-cone vector with f ^ = and Va ■ Pa = ^- If we work in a singular gauge, like 
the light-cone gauge, an additional gauge link at spatial infinity = +00) will have to be 
included in order to ensure the gauge invariance of the above definitions 14 1. 

Similarly, the quark distribution for the unpolarized hadron B can be defined as 

= ^-[q^'^'Hx,k^h,P^ + ...] , (9) 

where the gauge link is Cvt{g2', 00; ^) = exp [—ig2 dXvb ■ A{\v + ^)), with Vb a vector 
conjugate to Pb satisfying f ^ = and Vb ■ Pb = 1- Notice that in the above definition the 
gauge link is associated with the coupling g2. 

At the lowest order, in the one-gluon contribution case, the relevant hard factors in this 
model read: 

f^Sivers _ ghl gi+'2g22{s^ + U^) ^^^^ 
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FIG. 1: Two-gluon exchange contributions proportional to gf. 



where s, t, u are the partonic Mandelstam variables. The above results can be directly 
obtained from our previous result in 3] by replacing there C„ by 1, Ci and Cf2 by 92/91, 
and Cpi by 1. The expressions for the hard-scattering factors for other partonic channels 
may be found in a similar way. We note that the above result for H^gJ^^g^^, also agrees with 



that given in [8|. 

Before we move on to the two-gluon exchange contributions, we note that we have in 



fact also adapted our previous calculations of 



m 



l^see also 



20]) to this model case. As 



, LZj , we have started from a coUinear-factorized framework and considered the one-gluon 
contributions in the intermediate transverse momentum region Aqcd ^ 5'± ^ -Pi- We 
have found that indeed the factorized structure of Eqs. (I5l),(l6]) emerges, with perturbatively 
defined TMD distributions that depend on the charges 91 and 92 as given in Eqs. ([7]), ({9]), 
and with the hard-scattering factors in ffTOl) . We refrain from giving the details of these 
calculations in this paper, but stress that the one-loop TMD factorization for the dijet- 
imbalance found in [a, |7|] therefore extends to the model of 8|] we consider here. 

We now turn to the two-gluon exchange diagrams with the spectator of hadron A, and 
their contributions to the SSA and the spin-averaged cross section. They can be grouped 
into three parts: one proportional to 9f as shown in Fig. [1], where both gluons attach to the 
outgoing quark line with charge 91; one proportional to 9192 as shown in Fig. [21 where one 
gluon attaches to 91 quark line and another gluon to the incoming or outgoing 92 quark line; 
and one proportional to 92 as shown in Fig. [3], where both gluons attach to the 92 quark 
line. 

Let us first discuss the contributions from the diagrams in Fig. [H In the following calcu- 
lations, we only consider the dominant contributions from the exchange of gluons collinear 
to the polarized hadron A, which survive at the leading power of q±/P±- We can then utilize 
the eikonal approximation 

mm 

. The contribution from the diagram in Fig. [H^a) will 
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FIG. 2: Two-gluon exchange contributions proportional to gig2- 
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(i) 



FIG. 3: Two-gluon exchange contributions proportional to g2- 



then be proportional to 



9l 



1 



1 



(27r)4 (27r)4 (P4 - A;' - fci - ^2)2 - m2 + ie {k' + fci)^ - /x^ + (fc' + ki + ^2)^ - /^^ + 
1 1 / 1 \ / 1 



X 



fcf - A2 + kl-X^ + it \ -A;+ + ie 7 V -^1" - k'^ + ie 
where ki and A;2 are the momenta of the two exchanged gluons, and where we have introduced 
the masses m, /x, A for the scattering quarks and for the gluon, respectively. The last two 
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factors in (fTT!) come from the eikonal approximations for the two gluon attachments to the 
outgoing quark hne gi. A similar contribution is obtained from diagram (b), except for a 
difference in the eikonal propagators, for which we now have 

I \ f 1 



(12) 



V2 -r «e / \—k^ — k2+it^ 

Thus, the total contribution by these two diagrams will depend on the following expression 
for the eikonal propagators, 

1 \ / 1 \ f I \ ( 1 



A = gt 



-k{ + ie 



(13) 



and the rest of the expression is identical for these two diagrams. The imaginary part 
of the above expression contributes to the SSA, whereas the real part contributes to the 
unpolarized cross section. We can further simplify the above term as 

1 



A = gt 



+ tng^ 



kt 



(14) 



Similar calculations can be performed for the diagrams in Figs. [2] and [31 The contributions 
from diagrams 2 (a,b,d,e) will be proportional to the factor 

1 \ f 1 \ f 1 \ / 1 



B 



9i92 



+ 



—k^ + ie J \—k2 + ie 
1 \ / 1 



k^ + ie J \ —k2 + ie 
9i92 [4:{-i7T)^6{k+)6{k+)] +i27Tgig2 



—k^ + ie J \k2 + ie 



(15) 



where partial cancellations between the diagrams have occurred in the final result. The 
contributions from diagrams 3 (a,b,d,e,g,h) will depend on the factor 

1 \ / 1 \ / 1 \ / 1 



C 



9l 



+ 



+ 



—k'l + ie J \k2 + ie 
I \ f 1 



+ 



k{ + iej V^i + ^: 



le 



k{ + ie 



+ 



le 



1 



1 



-ki + ie J \ —kf — k2 + ie 
= gl [i{-t7ry5ikt)5{kt)] . 



k2 + ie J \k]^ + k2 + ie 
1 \ / 1 



+ 



le 



(16) 



Again, there are cancellations between the diagrams, which in this case leave us with only 
a real part contribution. 



The total contribution from the diagrams discussed above thus becomes 



A + B + C = gl 



1 



+ {-n,)H{k+)5{kt) + gi{gi + 2g2){iiT) 



6{kt) , 6{kt) 




(17) 



As we mentioned above, the first and third terms will contribute to the unpolarized cross 
section, whereas the second term will contribute to the SSA. Clearly, the first term can 
be factorized into the unpolarized TMD quark distribution, multiplied by the leading-order 
hard factor H^^,^^^, in Eq. ( fTOl) . The second term can also be factorized into the quark Sivers 
function for hadron A, multiplied by the leading-order hard factor H^^T^^^q, in Eq. ( fTOl) . The 
third term, however, cannot be factorized into the unpolarized quark distribution defined in 
Eq. ([7]), multiplied by the leading hard factor H^g,_^gg,. 

In order to account for this uncancelled contribution, we have to modify the gauge link 
definition in Eq. ([8]): 



which corresponds to what has been proposed in [2, y] . The first term takes into account 
the contributions due to final-state interactions of the outgoing quark with charge gi, the 
second one those of the outgoing quark with charge g2, and the third one the initial-state 
interactions of the incoming quark with charge g2- With this modification of the gauge link 
in the definition of the TMD quark distributions, the above results can all be reproduced by 
the two-gluon exchange diagram contributions, multiplied by the leading order hard factor. 
This is true for both the unpolarized and the spin-dependent (Sivers) case, even though for 
the latter, as seen from Eq. (ITTl) . the complications that require the redefinition of the gauge 
link are not yet visible at this order but should first occur at the next order of perturbation 
theory. 

A similar analysis can be performed for the diagrams in Fig. [l](c). Fig. [2]^c,f), and 
Fig. [3](c,f,i). All their contributions are reproduced by the above modified TMD quark 
distribution, multiplied by the leading-order hard factor. It is important to note that these 



^vaiduoo;^) C'^^{gi,g2;0 = V exp [ -igi 




(18) 
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diagrams do not contribute to the SSA. Moreover, their contributions, along with that in 
the last term of Eq. f|T7|) . will lead to an infrared- finite contribution to the differential un- 
polarized cross section. It will be interesting to investigate further how this will affect the 
dijet-correlation at hadron colliders. 

In the QED-like model of Ref. [8j], further study shows that we can simplify the gauge 
link in Eq. ([18]) to 

(0 = V exp (^-igi ^ dXva ■ A{i + Xva)^ 

xVexp(^-ig2j dXva-Ai^^ + Xva)^ ■ (19) 

We notice that the last factor will reduce to a pure phase when integrated over the transverse 
momentum, because in that case C,± is set to zero and there is no other dependence on ^. 
In this coUinear factorization approach is appropriate, and the gauge link will be 

defined in a standard way 2l|. In particular, the last term of Eq. (fT71) will be cancelled out 
by diagrams (c) and (c,f) of Figs. [2] and [31 respectively. 

4. ^Conclusions. In this paper, we have compared three recent theoretical studies 

m 

of the factorization properties of dijet-correlations in a polarized hadronic reaction 
when the two jets are produced nearly back2,to-back in azimuthal angle. We have shown 

inks derived 
st, we have 



[fl 



that the first-order perturbative results of [y, [3] are reproduced when the gauge 



m 



3| are expanded to that order. Within the abelian model considered in 
also calculated the two-gluon contributions to the initial-state and final-state interactions 
and verified that the results are consistent with the expected gauge-link structure. We 
therefore conclude that the results of the various studies are mutually consistent. They 
show that beyond the first order, the TMD quark distributions need to be modified in 
a process-dependent way, in order for them to correctly take into account coUinear gluon 
contributions. This is a departure from standard factorization, which on the other hand still 
does not exclude that the dijet observable obeys a more generalized factorization in terms 
of more complicated TMD correlators. Further studies will be needed here. 
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